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摘要 : 原子 钟 通过 稳定 的 频 标 来 测量 时 间 ， 但 频 标 输出 信号 会 受到 不 同 来 源 噪声 的 影响 ， 因 此 
原子 钟 间 相 位 偏差 可 以 看 作 服 从 突 律 SO ha f^ 的 连续 随机 过 程 。 该 随机 过 程 并 非 平稳 的 ， 一 般 
是 通过 高 阶 的 差分 使 得 序列 平稳 化 ， 常 用 的 时 域 频 率 稳定 度 评估 方法 便 是 采用 此 种 思路 。 血 律 噪 
声 实际 是 对 白 噪 声 过 程 进行 微分 和 积分 得 来 ， 不 同 肾 律 的 噪声 对 应 不 同 的 随机 微分 方程 ， 实 际 上 
窜 律 噪声 的 仿真 也 是 通过 白 噪声 基于 此 实现 的 。 首 先 介绍 了 不 同 虹 律 噪声 的 鉴定 方法 ， 并 结合 文 
献 给 出 常见 寡 律 噪声 对 应 的 微分 方程 ， 然 后 介绍 了 不 同 频率 稳定 度 评估 方法 之 间 的 关系 和 其 相 
应 的 传递 函数 ， 并 简单 总 结 了 其 置信 区 间 的 计算 方法 。 该 文 的 工作 有 助 于 构建 原子 钟 随机 模型 和 
频率 稳定 度 评估 方法 。 
X 键 词 : 原子 钟 ， 鹤 律 谱 噪 声 ， 等 效 自由 度 ; 置信 区 间 ; 无 偏 估计 
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1 引 8 


原子 钟 通过 稳定 的 频 标 来 测量 时 间 ， 最 终 输 出 的 钟 面 时 不 可 避免 地 受到 电子 元 器 件 运 行 
噪声 影响 ， 不 同 层 级 的 白 噪 声 随 着 时 间 积 累 体 现在 钟 面 时 的 功率 谱 密度 上 呈现 过 律 ^ ha f^ 
特性 ， 包 含 不 同 的 周期 特征 ”"。 自 相关 函数 (autocorrelation function, ACF) 和 功率 谱 密度 
(power spectral density, PSD) 包含 一 个 连续 随机 过 程 的 全 部 信息 ”， 但 很 难 从 中 提取 有 用 
的 定量 信息 。 阿 伦 方差 (Allan variance) ^ 595353477; 2: (Hadamard variance)” 可 以 定量 地 
描述 原子 钟 频率 稳定 度 ，Hutsell” 给 出 两 者 与 卡尔 曼 滤波 模型 过 程 噪声 系数 的 数学 关系 。 阿 
伦 方差 可 以 识别 不 同 虹 律 的 频率 白 噪声 、 频 率 内 变 噪声 和 频率 游 走 噪声 (对 于 频率 数据 ， 其 
RES) BIN 0, —1, -2)， 并 借 此 计算 噪声 水 平 系数 ”。 但 对 频率 游 走 闪 变 噪声 CREE —3) 
和 频率 奔跑 噪声 (WEA -4) 不 收敛 "”， 不 幸 的 是 鲍 原 子 钟 受 到 这 两 类 噪声 的 影响 ”。 与 阿 
伦 方差 不 同 ， 基 于 三 次 差分 的 哈达 玛 方差 对 这 两 类 低频 噪声 收敛 ”。 
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当 平滑 时 间 较 长 时 ， 哈 达 玛 方差 估 值 置信 度 会 降低 ， 有 学 者 提出 一 种 总 方差 的 估 值 方 
法 ， 该 方法 通过 特定 的 映射 来 延 拓 数 据 以 提高 置信 和 度 ” "。 参 考 重 又 阿 伦 方差 ”类 似 地 可 以 
定义 重 肥 哈达 玛 方差 "， 由 于 重 肥 哈达 玛 方差 最 大 限度 地 利用 了 样本 的 所 有 可 能 的 三 次 差 ， 
因而 最 大 限度 地 利用 了 当前 样本 ， 所 以 其 估 值 的 置信 度 可 以 得 到 提升 。 与 哈达 玛 方差 的 估 值 
不 同 ， 重 闭 哈 达 玛 方差 采用 了 彼此 存在 相关 性 的 样本 进行 估 值 ， 也 因此 其 估 值 不 服从 整数 
自由 度 的 卡 方 分 布 ”， 而 是 经 验 性 地 近似 服从 非 整 数 等 效 自由 度 v 的 卡 方 分 布 “。 对 于 方 
差 的 估 值 ， 当 显著 性 水 平 相同 时 ， 等 效 自由 度 越 大 估 值 置信 区 间 越 小 ， 估 值 越 准确 。 随 着 北 
斗 三 号 全 球 系统 的 建立 ， 对 高 性 能 星 载 原子 钟 钟 差 的 计算 、 预 报 和 评估 工作 得 以 开展 77, 
而 开展 此 类 工作 的 前 提 便 是 构建 这 些 原子 钟 的 基本 数学 模型 。 

本 文 第 2 章 首 先 介 绍 了 最 广 为 接 受 的 原子 钟 模型 ， 该 模型 可 以 说 是 最 简单 有 效 的 原子 
钟 模型 ， 其 同时 包含 了 确定 过 程 和 随机 过 程 部 分 ， 其 随机 过 程 部 分 即 为 非 平 稳 的 宕 律 谱 模 
型 。 对 于 原子 钟 频 率 而 言 ， 一 般 认 为 其 受到 5 种 整数 过 律 噪声 (2 ~ -2) AMM", Ju 
续 对 锦 原 子 钟 的 研究 增加 了 两 种 寡 律 噪声 (—3, -4) ”， 本 文 第 2 章 还 介绍 了 不 同 窜 律 噪声 
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2 律 为 —2 的 随机 过 程 ， 由 于 其 为 独立 增 量 过 程 ， 对 其 微分 便 是 白 噪声 过 程 。 本 文 第 3 章 主 要 
c 介绍 了 窜 律 噪声 的 仿真 以 及 其 对 应 的 随机 微分 方程 ， 至 此 主要 为 时 域 分 析 。 原 子 钟 频率 稳定 


度 的 评估 便 是 对 容 律 噪声 能 量 的 评 佑 ， 阿 伦 方差 与 哈达 玛 方差 在 本 质 上 是 一 致 的 ， 即 在 时 
域 上 通过 差分 使 得 序列 平稳 化 。 本 文 以 概述 的 角度 在 第 4 章 统一 介绍 了 两 种 频率 稳定 度 评 
估 方 法 及 其 关系 ， 同 时 整理 了 它们 在 频 域 的 传递 函数 。 对 于 一 个 高 斯 过 程 ( 白 噪声 过 程 ) 而 
言 ， 最 重要 的 统计 量 便 是 期 望 和 方差 国 数 ， 其 期 望 常 被 认为 是 零 ， 其 方差 主要 反映 了 噪声 能 
量 的 强 弱 。 众 所 周知 ， 一 般 对 方差 估计 的 统计 量 服从 卡 方 分 布 ， 那 么 对 于 阿 伦 或 者 哈达 玛 方 
差 而 言 ， 其 置信 区 间 同 样 通过 卡 方 分 布 确定 。 第 5 章 主要 介绍 等 效 自由 度 和 频率 稳定 度 评 
佑 方法 置信 区 间 的 确定 。 


2 原子 钟 模型 和 窜 律 噪声 的 鉴别 


原子 钟 借助 稳定 的 频 标 测量 时 间 ， 通 过 比 对 原子 钟 间 的 相位 差 ， 根 据 频 标的 基准 频率 得 
到 物理 意义 上 的 时 间 差 。 但 原子 钟 频 标 很 容易 受到 不 同 来 源 的 噪声 影响 ， 温 度 和 辐射 同时 
会 影响 原子 钟 频 标的 稳定 度 ”。 一 种 被 广 为 接 受 且 相对 简单 的 原子 钟 频 标 电子 振荡 器 模型 
为 : 


V(t) = (Vo + &(t)) sin[27 + O(t)] , (1) 


EP, Vy Al vo ASAT AR, n 是 圆周 率 ，t 是 时 间 ; e(t) 是 施加 在 振幅 上 的 噪声 ， 
对 于 精确 的 原子 钟 振荡 器 而 言 ， 该 噪声 可 以 忽略 ; e(t) 是 一 个 非 平 稳 的 连续 窜 律 随机 过 程 
Hi [ó(t)| « 2r. EET POK 2rxo， 即 该 噪声 对 频率 的 影响 远 低 于 原子 钟 基准 频 
率 ， 更 容易 反映 在 经 过 一 定时 间 积 累 后 的 相位 和 频率 比 对 序列 上 。 

一 般 一 个 原子 钟 的 频 标 被 认为 是 一 个 十 分 稳定 的 频率 源 ， 即 频 标的 输出 频率 保持 在 基准 
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频率 vo 附近 很 小 的 一 个 区 间 内 。 原 子 钟 对 时 间 的 测量 通过 相位 比 对 实现 ， 当 相位 由 如 .2rmm 
变化 到 如 . 27g + At 2 = (fo + At) 27g 时 ， 时 间 便 增加 了 At 秒 。 但 一 个 频 标的 实际 
频率 必然 与 其 标 称 频率 之 间 存在 一 个 系统 差 ， 初 始 相位 也 往往 不 可 能 是 零 ， 同 时 还 可 能 伴 
随 着 频率 老化 和 频率 漂移 现象 ， 因 此 原子 钟 测量 的 时 间 偏差 会 受到 初始 相位 偏差 zo、 初 始 
频率 偏差 yy 和 线性 频率 漂移 D 等 确定 部 分 的 影响 。 同 时 原子 钟 的 相位 还 受到 随机 部 分 (t) 
的 影响 。 那 么 ， 瞬 时 时 间 (相位 ) 偏差 z(t) GE UN: 


Deystem Eu p(t) 
27t/9 


FEA, Deystem 代表 由 初始 相位 偏差 zo、 初 始 频率 偏差 yo。 和 线性 频率 漂移 D 等 确定 过 程 导 
致 的 相位 偏差 ，9(t) 是 由 随机 过 程 导致 的 相位 偏差 ， 其 功率 谱 是 Self) 20, RP Salf) 
表示 功率 谱 密度 。 

自然 环境 中 原子 钟 振荡 器 的 输出 频率 不 可 能 是 恒定 的 : (1) 频 标 输出 频率 本 身 存 在 白 品 
Pi. (2) 受 环境 因素 和 频 标 老 化 影响 ， 频 标 输出 频率 会 出 现 包 含 确定 过 程 的 增 量 过 程 ， 该 增 
量 过 程 也 会 受到 白 品 声 污染 并 使 得 频 标 输出 频率 包含 频率 游 走 噪 声 。 因 此 ， 最 初 的 原子 钟 随 
机 模型 仅 为 暴 律 2 到 -2 的 噪声 受 加 ， 其 中 肾 律 为 2 的 噪声 是 自然 存在 于 相位 比 对 时 的 白 噪 
E, TEE 0 的 噪声 为 频 标 输出 频率 本 身 存在 白 噪声 ， 需 律 为 -2 的 噪声 为 上 述 的 频率 游 走 
噪声 ， 早 律 为 1 和 一 1 的 噪声 被 称 为 内 烁 噪声 ， 其 由 电子 元 器 件 以 及 电子 环 路 非 理 想 因 素 
致 ”。 而 原子 钟 间 的 相位 偏差 z(t) 实际 上 是 原子 钟 频 率 偏差 的 积分 过 程 ， 因 此 有 : 


a(t) = 


1 
= £o + yot + 5 Dt Fest) , (2) 


m 


emn 


B eu eno s hu. (3) 


a=—4 


其 中 ，S,(f) 是 oz 的 的 功率 谱 密 度 (PSD), S,(f) 是 wy 人 的 功率 谱 密度 ， 且 有 plt) = 
g(t); 对 于 频率 f 0x f x fo KU fy 为 上 限 截止 频率 ;ea IRERE EHG ha 
ERAO a 的 噪声 的 强度 因子 。a = 2, 1,0,… ,一 4 分 别 对 应 噪声 名 为 相位 白 噪声 (white 
phase modulation noise, WPM)、 相 位 闪烁 噪声 (flicker phase modulation noise, FPM), 
频率 白 噪 声 (white frequency modulation noise, WFM), MiX NARRE (flicker frequency 
modulation noise, FFM)、 频 率 随 机 游 走 噪声 (random walk frequency modulation noise, 
RWFM)、 频 率 内 烁 游 走 噪声 (flicker walk frequency modulation noise, FWFM) 和 频率 随机 
奔跑 噪声 (random run frequency modulation noise, RRFM)。 前 五 种 噪声 便 是 上 文 提 到 的 
最 初 原子 钟 随机 模型 里 所 建 模 的 五 种 噪声 ， 后 两 种 噪声 是 因为 后 续 对 鲍 原 子 钟 研究 而 增加 
的 奔跑 闪烁 噪声 和 奔跑 噪声 。 其 中 奔跑 噪声 可 以 理解 为 “随机 游 走 的 游 走 噪声 ?”， 即 频 标 输 
出 频率 的 增 量 过 程 也 受到 了 白 噪声 的 污染 。 式 (3) 便 是 有 名 的 原子 钟 噪声 窜 律 模型 "。 

对 于 如 何 确 定 不 同 窜 律 噪声 的 类 型 。 最 为 直接 的 想法 就 是 估计 噪声 功率 谱 ， 从 频 域 上 直 
接 分 析 。 但 此 方法 并 不 容易 实现 ， 原 因 主 要 为 两 个 : (1) 实际 的 噪声 往往 是 多 窜 律 噪声 的 舍 
加 ; (2) 功率 谱 的 准确 估计 本 身 就 是 难题 。 一 般 采 用 时 域 分 析 的 方法 辨别 噪声 ， 这 里 介绍 一 
种 基于 时 域 的 自 相关 函数 的 噪声 鉴别 方法 ，Lag 1 偏差 函数 法 ”。 
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Lag k 自 相 关 函 数 (ACF) 的 定义 为 : 


py = Ele Wee a 


HP, 2, 是 任 一 时 间 序 列 ， 可 以 是 瞬时 相位 偏差 z( 或 瞬时 频率 偏差 y(t) = z(t) 的 离散 采 
样 序列 ; n 是 其 均值 ，o? BHAA, El] 代表 数学 期 望 。 其 为 交错 位 后 的 序列 自 相 关 函 
数值 。 当 大 = 工时 ， 便 是 Lag] 自 相 关 函 数 ， 定 义 统计 量 rk 作为 其 估计 值 : 


k > (21 — Z) (Zt+k — 2) 
Tk = ; (5) 


N DE — zy 


其 中 , z 是 该 时 间 序 列 的 均值 ，N 是 该 时 间 序 列 的 长 度 。 

TIEN 2 的 噪声 实际 上 是 白 品 声 过 程 经 过 一 次 微分 得 到 的 随机 过 程 ， 对 于 该 过 程 均匀 
采样 的 时 间 序 列 而 言 ， 相 邻 的 两 个 随机 变量 间 存 在 线性 相关 性 ， 其 Lag 1 自 相 关 函 数值 为 
-1/2， 该 值 可 由 维 纳 - 辛 钦定 理 推导 出 。 而 过 律 为 1 的 噪声 同样 是 白 噪 声 过 程 经 过 更 复杂 的 
微分 (分 数 微分 ) 得 到 ， 其 Lag 1 自 相关 函数 值 为 —1/3. VUE EE S 0 的 噪声 实际 就 是 白 噪 
声 ， 对 该 过 程 均匀 采样 的 时 间 序 列 相 邻 的 两 个 随机 变量 线性 无 关 ， 故 其 Lag 1 自 相 关 函 数值 
为 零 。 图 1 展示 了 仿真 的 占 律 为 2, 1 和 0 的 随机 过 程 的 20 组 样本 函数 (噪声 的 强度 因子 为 
1， 时 间 采 样 间隔 为 1 s) 的 自 相 关 函 数 图 ， 仿 真 的 Lag 1 自 相 关 函 数值 与 理论 分 析 一 致 。 


1 HEA 2, 1 和 0 的 随机 过 程 的 20 组 样本 函数 的 自 相关 函数 图 


E a 是 频率 偏差 序列 ，p1 的 估 值 71 接近 一 1/2 便 是 WPM、 接 近 一 1/3 便 是 FPM、 接 
近 0 便 是 WFM， 这 三 种 噪声 对 于 频率 偏差 序列 而 言 都 是 平稳 的 ， 可 以 直接 通过 Lag 1 自 相 
关 函 数 确定 其 噪声 类 型 。 对 于 非 平稳 的 时 间 序 列 ， 需 要 首先 通过 差分 使 其 平稳 化 进而 利用 
Lag 1 自 相 关 函 数 确 定 其 噪声 类 型 。 令 5 = r/r) 一般 的 若 5 < 0.25， 则 可 以 认为 时 
间 序 列 z, 平稳， 并 有 Lag 1 相关 函数 p = 0/(1 一 5)。 这 样 ， 对 于 一 个 非 平稳 的 时 间 序 列 
a: 首先 对 其 进行 多 次 向 后 一 阶 差 分 并 计算 6 直到 6 < 0.25， 此 时 时 间 序 列 平稳 。 注 意 到 每 
差分 一 次 得 到 新 的 序列 ， 对 于 新 序列 来 讲 原 本 噪声 的 容 律 +2， 此 时 记录 差分 次 数 d 用 以 还 
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原 原 本 噪 声 宕 律 系数 。 经 过 差分 后 的 序列 实际 是 把 寡 律 为 -1 ~ -4 的 噪声 转化 为 寡 律 为 2 
~ 0 的 噪声 ， 此 时 利用 Lag 1 ACF 函数 确定 具体 的 噪声 类 型 并 结合 差分 次 数 d 得 到 原本 品 
ERIAM S p= -round(26) — 2d (其 中 round(*) 为 取 * 最 近 的 整数 )， 可 借助 p 的 值 
Hf iE BRE EHE w。 对 于 频率 数据 而 言 a = p， 对 于 相位 数据 而 言 a= p42 

常用 的 宕 律 噪声 鉴定 方法 还 有 Bl 偏差 函数 法 ，” 等。 对 于 Lag 1 偏差 函数 法 ， 文 献 [30] 
做 了 进一步 拓展 ， 并 给 出 了 基于 重 又 哈达 玛 方 差 的 Lag 1 自 相 关 函 数 守 律 噪声 确定 算法 。 


3 ”和 窜 律 噪声 的 仿真 及 其 随机 微分 方程 


白 噪声 过 程 的 功率 谱 密 度 S, (u) 在 整个 频谱 上 为 常数 ， 对 和 白 噪 声 过 程 w(t) 进行 积分 和 
微分 就 能 生成 频谱 为 暴 律 的 有 色 品 声 过 程 。 帘 律 为 2, 0, -2 和 -4 的 有 色 噪 声 可 通过 白 噪声 
过 程 直接 微分 或 者 积分 得 到 : 


(po =w 
dy-2 ; 6 
m (6) 
d^. 4 
dt? 
Hh. o; ERIEN d PR ERES, w 为 白 噪声 过 程 。 式 (6) 给 出 窜 律 为 2, 0, —2 和 —4 等 有 
色 噪 声 的 随机 微分 方程 ， 借 助 该 方程 可 以 由 和 白 噪 声 过 程 仿真 出 上 述 有 色 噪 声 。 实 际 上 最 为 ) 
泛 的 基于 卡尔 曼 滤波 实现 的 时 间 尺 度 的 状态 方程 便 建立 在 此 基础 上 。 式 (7) 为 标准 三 变量 卡 
尔 曼 滤波 时 间 尺 度 用 以 描述 系统 状态 的 微分 方程 组 : 


= wW 


— 


dx, 

II ee +w 

di T2 1 

dx 

m = T3 T W2 ’ (7) 
Oy 

dt ^ 5 


HP, 21, to 和 zs 分别 为 相位 状态 量 、 频 率 状态 量 和 频率 漂移 状态 量 ，wi, w 和 ws 为 三 个 
独立 的 白 噪声 过 程 ， 其 分 别 生成 了 FFM, RWFM Ñ RRFM””. 

SAM TEA 1, —1 和 -3 的 有 色 噪 声 ， 仿 真一 直 是 个 难题 ， 主 要 的 仿真 方法 有 三 
FH: (1) 借助 ARMA 模型 仿真 ; (2) 分 数 微 积分 仿真 ; (3) 基于 小 波 分 析 仿真 。 三 种 仿真 方 
法 都 需要 采用 不 同 的 有 限 过 程 逼 近 寡 律 为 1, —1 和 —3 的 有 色 噪 声 汪 。 文 献 [27] 对 这 类 噪声 
产生 的 机 理 进 行 了 深入 探 帘 ， 并 找到 了 产生 此 类 噪声 背后 的 物理 原因 。 本 文 仅 结合 文献 介绍 
MEN -1 的 有 色 噪 声 定 义 ， 以 及 采用 四 个 马尔 可 夫 过 程 逼近 的 噪声 生成 方法 。 考 虑 到 过 律 
为 —2 的 噪声 就 是 布朗 运动 ， 借 助 修正 的 分 数 布朗 运动 : 

1 


Bj) = aay || (emn mieten f -rular ` (8) 
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定义 任意 a ARR. KEREI —1 的 噪声 在 数学 上 便 是 采用 分 数 微 分 定义 的 ， 
即 取 公 式 (8) 中 的 a = 1。 然 而 仅 借 助 公式 (8) 不 可 能 实现 分 数 布朗 运动 仿真 ， 本 文 介绍 一 
种 基于 多 个 独立 马尔 科 夫 过 程 又 加 的 究 律 为 —1 的 噪声 仿真 方法 。 

一 个 马尔 可 夫 频 率 过 程 可 由 高 斯 过 程 借助 微分 方程 化 ， 表 示 为 : 


— + hom =w , (9) 


其 中 ，ww 是 马尔 可 夫 频 率 过 程 ，R 是 常数 ，w 是 白 噪声 过 程 。 选 择 不 同 的 常数 R 生成 四 个 
独立 的 马 氏 过 程 并 琶 加 ， 便 可 近似 地 生成 一 个 寡 律 为 -1 的 有 色 品 声 过 程 ”。 文 献 [34] 同 
时 完成 了 包含 FFM 噪声 的 卡尔 曼 时 间 尺 度 的 构造 ， 其 对 FFM 噪声 的 处 理 办 法 便 是 采用 四 
个 独立 的 马 氏 过 程 近似 ，R 作为 状态 参数 被 设计 在 滤波 器 中 用 以 拟 合并 逼近 原子 钟 里 实际 
FER FFM. 

图 2 展示 了 白 噪声 过 程 (REIO Ass) EU -—1). AB) Efe —2) 
和 随机 奔跑 运动 (EEA —4) 4 种 随机 过 程 的 20 组 仿真 样本 时 间 序 列 图 。 借 助 该 图 可 以 直 
观 地 了 解 不 同 暴 律 噪声 在 时 域 里 的 特征 。 


10}: i a=-1 
m | 
E 
| « &§ i 
pra OR E Qo erna TITAN 
1000 2000 3000 4000 0 1000 2000 3000 4000 
时 间 /s 时 间 /s 
a) b) 


仿真 噪声 值 


78008 1000 2000 3000 4000 -50 1000 2000 3000 4000 
时 间 /s 时 间 /s 
c) d) 


HE: a) 白 噪 声 过 程 ; b) 分 数 布朗 运动 ，c) 布朗 运动 ; d) 随机 奔跑 运动 。 其 中 仿真 的 噪声 强度 因子 为 1， 样 本 
函数 时 间 采 样 间隔 为 1 so 


2 ”连续 随机 过 程 的 20 组 样本 函数 时 间 序 列 图 


原子 钟 系统 模型 和 随机 模型 已 梳理 清楚 ， 帘 律 噪声 的 鉴别 与 仿真 已 完整 介绍 。 那 么 此 后 
借助 常用 的 频率 稳定 度 评估 方法 可 以 进一步 从 多 种 混合 的 帘 律 噪声 中 确定 相应 噪声 的 容 指 
数 和 能 量 系数 。 
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4 原子 钟 频率 稳定 度 评估 方法 的 统一 形式 和 传递 函数 


常用 的 原子 钟 频 率 稳定 度 评估 方法 有 阿 伦 方差 和 哈达 玛 方差 ， 重 又 阿 伦 方差 和 重 又 哈 
达 玛 方差 是 为 解决 数据 量 少时 两 种 方差 置信 度 低 的 问题 ， 改 进 阿 伦 方差 和 改进 哈达 玛 方差 
是 因为 两 种 方差 无 法 区 分 寡 律 为 2 和 1 的 噪声 ， 而 阿 伦 总 方差 和 哈达 玛 总 方差 是 为 解决 平 
滑 时 间 较 长 时 置信 和 度 低 的 问题 所 提出 ”。 这 些 频 率 稳定 度 评估 方法 都 是 基于 阿 伦 方差 和 哈 
达 玛 方差 ， 本 章 主 要 讨论 两 种 方差 的 定义 、 关 系 及 传递 函数 ， 并 结合 文献 介绍 一 种 统一 的 形 
式 用 来 描述 两 种 方差。 
阿 伦 方差 为 分 析 频 率 稳 定 度 最 早 被 提出 ， 其 本 质 仍 是 一 种 方差 ， 由 于 振荡 器 频率 序列 受 
到 随机 游 走 噪声 影响 并 非 平稳 序列 ， 因 此 需要 通过 差分 使 得 序列 平稳 化 。 频 率 稳定 度 实际 评 
估 的 是 一 个 振荡 器 输出 频率 的 随机 误差 ， 与 方差 分 析 概 念 有 些 类 似 但 内 涵 完 全 不 同 ， 这 是 由 
基本 的 随机 模型 不 同 所 决定 的 。 频 率 稳定 度 评估 的 随机 模型 是 一 种 不 平稳 的 随机 过 程 ， 是 
多 种 不 平稳 过 程 的 三 加 ， 我 们 需要 对 这 个 多 种 噪声 玛 加 的 不 平稳 过 程 进行 整体 的 误差 分 析 。 


Tr 


T 而 测量 平 差 中 的 误差 分 析 其 实 是 对 平稳 的 白 噪声 的 误差 分 析 。 当 然 两 者 也 存在 一 定 的 联系 ， 
© 当 我 们 对 WFM, RWFM 和 RREM 进行 建 模 ， 引 入 状态 空间 ， 此 时 在 状态 空间 中 的 误差 分 
析 便 转化 为 了 测量 平 差 中 的 误差 分 析 。 这 种 转化 实际 上 是 把 “基于 平滑 时 间 变 化 的 多 品 声 车 


加 随机 误差 分 析 ” 转 化 为 了 “状态 空间 中 状态 量 的 误差 分 析 ” 但 事实 上 ， 对 某 一 种 究 律 品 
声 的 解析 式 的 独立 分 析 本 质 上 是 很 难 实现 的 。 对 原始 数据 进行 差分 后 ， 数 据 的 物理 意义 已 
经 改变 ， 噪 声 的 功率 谱 也 发 生 了 变化 ， 但 频率 稳定 度 所 分 析 的 目标 是 一 个 频 标 输出 的 随机 
误差 ， 是 对 晕 律 谱 随 机 模型 的 整体 分 析 ， 差 分 后 噪声 功率 谱 确实 发 生 了 一 些 变化 ， 由 Sf) 
变 成 了 (2xf)?. 5S(f)， 但 这 并 不 影响 对 频 标 输出 的 随机 误差 的 整体 分 析 。 阿 伦 方差 便 是 为 分 
析 频 率 稳定 度 提出 的 ， 对 于 频率 数据 而 言 ， 由 于 只 有 作为 主要 噪声 源 之 一 的 WEM 的 方差 
是 不 随时 间 变 化 的 ， 因 此 阿 伦 方差 被 设计 为 其 期 望 与 WEM 的 标准 方差 的 一 致 。 阿 伦 方差 


c2(r) 的 定义 是 : 
sg) = (Bie P), ot f yoa, (10) 


T 


= KB, y(t) 是 瞬时 频率 偏差 ，r 是 平滑 时 间 。 由 于 原文 是 定义 在 连续 过 程 上 的 ， 需 要 在 
o t = [0, 十 co] 上 积分 ， 故 () 表示 无 穷 时 间 上 的 平均 ， 可 以 理解 为 数学 期 望 ""。 若 振荡 器 频率 
受到 随机 奔跑 噪声 影响 ， 仅 通过 一 次 差分 无 法 使 得 序列 平稳 化 ， 很 容易 想到 进行 高 阶 差分 使 
序列 平稳 化 。 实 际 上 哈达 玛 方差 是 基于 哈达 玛 变换 构造 的 ”， 文 献 咎 将 其 改造 使 其 更 加 适 
合 时 域 分 析 。 实 际 上 根据 文献 [4] 的 定义 ， 对 频率 数据 进行 二 阶 及 二 阶 以 上 (相位 数据 三 阶 
及 三 阶 以 上 ) 差分 构造 的 方差 均 称 为 哈达 玛 方差 ， 不 过 对 于 原子 钟 频率 稳定 度 分 析 来 讲 ， 对 
频率 数据 二 阶 以 上 的 差分 并 无 必要 ， 我 们 通常 讲 的 哈达 玛 方差 便 是 对 频率 数据 进行 二 次 差 
分 的 方差 。 哈 达 玛 方差 定义 为 : 


= — 4 +9 2 7 1 th+r 
oi (2,7) = (Gea Phe 997), a on f" var (11) 
tk 
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公式 (11) 可 以 拓展 为 N 次 差分 哈达 玛 方差 : 
E 7 2 p T 

oh (yr) = (Sein ect Rea ty, tf y(t)dt , (12) 
当 N = 1 时 ， 为 阿 伦 方差 N > 2 时 ， 是 哈达 玛 方差 。 需 要 注意 的 是 ， 式 (12) 的 定义 与 
文献 [4] 最 初 的 定义 不 太一 致 ， 首 先 差分 方式 略 有 不 同 ， 其 次 原本 的 定义 并 未 除 以 分 母 的 部 
分 。 考 虑 到 差分 运算 会 放大 噪声 ， 且 为 了 保证 该 方差 的 期 望 与 WFM 标准 方差 的 期 望 一 致 ， 
式 (12) 分 母 是 必要 的 。 式 (11) 所 示 的 现在 最 常用 的 哈达 玛 方差 是 经 过 了 Riley 和 Howe 改 
进 的 ， 他 们 同时 给 出 了 上 述 二 阶 差分 哈达 玛 方差 的 传递 函数 ” 。 另 外 ， 式 (12) 的 拓展 是 作 
者 根据 文献 归纳 推导 出 的 ， 但 遗憾 的 是 本 文 未 能 给 出 NN > 2 时 的 传递 函数 的 统一 形式 。 
通过 阿 伦 方差 和 哈达 玛 方差 的 定义 ， 不 难看 出 两 种 方差 的 计算 都 需要 经 过 差分 和 平滑 
两 种 算 子 的 处 理 ， 而 后 续 对 两 种 方差 的 改进 也 无 非 是 多 增加 平滑 。 众 所 周知 ， 平 滑 算 子 可 以 
看 作 一 种 低 通 滤波 器 ， 而 差分 算 子 可 以 看 作 一 种 高 通 滤波 器 ， 两 者 组 合 可 以 构成 一 个 带 通 
滤波 器 ， 并 且 随 着 平滑 时 长 + 的 增 大 ， 该 带 通 滤 波 器 的 最 大 共振 频率 随 之 降低 。 而 方差 本 
身 便 是 对 噪声 能 量 强 弱 的 一 种 评估 ， 因 此 通过 平滑 时 长 r 的 增加 ， 不 同 窜 律 噪声 将 逐次 变 
为 主导 噪声 类 型 ， 进 一 步 可 以 顺 次 分 析 WEM, FFM, RWFM, FWFM 和 RREM 等 噪声 的 
能 量 强 弱 。 式 (13) 和 (14) 分 别 给 出 了 阿 伦 方差 和 哈达 玛 方差 的 传递 函数 ”， 其 中 n (f) 和 
Hyl f) 分 别 是 阿 伦 方差 和 哈达 玛 方差 的 传递 函数 ，7 是 平滑 时 长 。 


mop = [Tem]. (13) 
ssp =2 [EE (14) 


5 方差 的 分 布 和 等 效 自由 度 的 确定 


阿 伦 方差 和 哈达 玛 方差 是 频率 稳定 度 评 佑 的 主要 工具 ， 而 这 两 种 工具 的 本 质 仍 是 一 种 
对 方差 的 估计 。 方 差 的 估 值 一 般 是 服从 卡 方 分 布 ， 那 么 对 方差 估 值 的 误差 分 析 就 应 当 基 于 卡 
方 分 布 而 非 正 态 分 布 。 本 章 主要 介绍 基于 卡 方 分 布 的 方差 估 值 误差 棒 的 计算 。 

对 于 一 个 平稳 过 程 方差 o? 的 无 偏 估 计 V (o? = Ve) FAC FARE (equivalent degrees 
of freedom, EDF) 73: 


2(Vg)? 

Var 
随机 变量 (v/o?) V 近似 地 服从 自由 度 为 v 的 卡 方 分 布 ， 那 么 给 定 的 显著 性 水 平 p 可 以 确 
定 方差 估 值 o? 的 置信 区 间 为 vV/zs < 0? < vV/o, (£1 < x2)。 文献 [14] 给 出 了 Ve 和 Voar 


v = Vgpr = 


(15) 


Ve = 0° =s,(0) , (16) 
2 M 
Var = M2 `> s?[(na BH ni)r] ) (17) 
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Kp, +r 为 平滑 时 长 。 自 相关 函数 s:( 由 广义 自 相 关 函 数 (generalized autocovariance, 
GACV) su(b) 表示 : 


s.(t) 一 (AcA. SALA. cst) ’ (18) 


其 中 ，A 为 一 阶 的 向 后 差分 算 子 ，qd = 2 时 为 阿 伦 方差 ， 4 = 3 时 为 哈达 玛 方差 。 对 于 不 同 
PERRE, s(t) 不 同 ， 表 128618 T 381873 2, 1, 0, 一 1, 一 2, 一 3, 一 4 的 广义 自 相 关 函 数 。 


a 2 1 0 1 2 3 4 
-lt mt Jj? -tnt -It t$mjt.— [u^ 


So (t, a 


wa 


按照 本 文 第 2 章 介绍 的 方法 ， 可 以 确定 平滑 时 长 r NS ER RAR, 23541 
并 借助 式 (16) 和 式 (17) 可 以 计算 出 方差 估计 的 等 效 自由 度 ”>， 那 么 该 方差 的 分 布 函数 为 : 


2 
V:s 
X% = 一 一 , (19) 


KH, s 是 采样 方差 ，o? 为 真实 方差 ，v 为 等 效 自由 度 。 给 定 显著 性 水 平 p， 双 边 置信 区 
间 为 : 


2 v 2 2 v 
(ay) 7^ BID 
ZTE SB ESL FARSI. 
基于 本 章 介绍 的 算法 ， 可 以 计算 阿 伦 方差 和 哈达 玛 方差 估 值 的 误差 棒 。 


(20) 


6 结 ie 


本 文 对 原子 钟 频率 稳定 度 分 析 方 法 的 数学 原理 进行 了 简单 的 介绍 。 上 有 具体 包 括 : 原子 钟 模 
型 的 建立 及 过 律 谱 噪声 的 识别 ， 才 律 噪声 的 仿真 及 其 随机 微分 方程 ， 常 用 频率 稳定 度 评估 方 
法 的 统一 描述 形式 及 其 传递 函数 ， 常 用 频率 稳定 度 评 佑 方法 的 置信 区 间 计 算 。 本 文 重点 介 
绍 了 频率 稳定 度 评估 算法 的 来 历 和 其 背后 的 数学 原理 。 通 过 本 文 的 综述 ， 可 以 对 原子 钟 随 
机 模型 和 频率 稳定 度 评估 算法 的 构建 得 到 系统 且 相 对 深入 的 认识 : (1) 究 律 噪声 是 什么 ，(2) 
最 初 的 原子 钟 随机 模型 为 何 仅 涉 及 五 种 肾 律 噪声 ，(3) 窜 律 噪声 如 何 鉴别 ，(4) TERR PS 
何 通 过 白 噪声 生成 ，(5) 不 同 频 率 稳定 度 评估 方法 的 内 在 联系 ; (6) 阿 伦 方差 和 哈达 玛 方差 
估 值 的 置信 区 间 如 何 计算 。 值 得 注意 的 是 ， 本 文 还 同时 给 出 一 种 包含 阿 伦 方差 、 哈 达 玛 方 
差 和 更 高 阶 差分 方差 统一 形式 。 基 于 本 文 介绍 的 相关 理论 ， 对 于 一 个 原子 钟 频率 稳定 度 的 
完整 评估 应 当 包 括 两 部 分 : (1) 确定 一 个 原子 钟 振荡 器 输出 频率 主要 受到 何 种 随机 噪声 的 影 
Wi; (2) 确定 影响 原子 钟 频率 稳定 度 的 几 种 主要 噪声 类 型 的 能 量 强 弱 。 基 于 本 文 介绍 的 相关 
理论 ， 对 于 相位 数据 或 者 频率 数据 ， 完 整 的 评估 过 程 应 当 为 : (1) 计算 在 不 同 平 滑 时 间 下 的 
阿 伦 方差 、 哈 达 玛 方差 和 改进 阿 伦 方差 ， 并 计算 方差 的 置信 区 间 ; (2) 计算 Lag 1 自 相 关 函 
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数 或 者 B1 偏差 函数 ， 确 定 不 同 平滑 时 长 下 的 主要 噪声 类 型 ，(3) 基于 此 完成 主要 噪声 类 型 
的 分 析 和 噪声 能 量 强度 分 析 。 
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A. Compilation of the Theoretical Fundamentals of Atomic 
Clock Frequency Stability Evaluation Methods 


MA Yue-xin!?, TANG Cheng-pan!?, HU Xiao-gong! 


(1. Shanghai Astronomical Observatory, Chinese Academy of Science, Shanghai 200030, China; 
2. University of Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: The atomic clocks measure time by the stable frequency standard, and the 
output signal of the oscillators is easily affected by noise from different sources which makes 
the phase deviation between two atomic clocks obeying a continuous power-law random 
process. Such stochastic process is not stationary and a typical approach by utilizing a 
high-order backward differencing operator can make the process become stationary, and the 
commonly used method for assessing the frequency stability in the time domain is based 
on this idea. This paper focuses on the mathematical fundamentals behind the method for 
assessing the frequency stability of atomic clocks. The power-law noise is actually derived 
from the calculus of white noise processes and the different power-law noises correspond 
to the different stochastic differential equations. In fact, the simulation of the power-law 
noise is also achieved by some differentiations and integrations operating on the white noise 
process. This paper first introduces the identification method of different power-law noises 
and gives the differential equations corresponding to the power-law noises in conjunction 
with some available literatures. Then, the relationships between the different frequency 
stability assessment methods and their corresponding transfer functions are presented, and 
the calculation method of the confidence intervals is briefly summarized. This paper provides 
some insight into the construction of the stochastic model for atomic clocks and methods for 


assessing frequency stability in a relatively systematic basis. 


Key words: atomic clocks; power-law spectral noise; equivalent degrees of freedom; confi- 


dence intervals; unbiased estimates 


